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Progress on Perfect Lattice Actions for QCD 

W. Bietenholz a , R. Brower a b , S. Chandrasekharan a and U.-J. Wiese a * 
a CTP-LNS, Massachusetts Institute of Technology, Cambridge MA 02139, USA 
b Dept. of Physics, Boston Univ., 590 Commonwealth Ave., Boston MA 02215, USA 

We describe a number of aspects in our attempt to construct an approximately perfect lattice action for QCD. 
Free quarks are made optimally local on the whole renormalized trajectory and their couplings are then truncated 
by imposing 3-periodicity. The spectra of these short ranged fermions are excellent approximations to continuum 
spectra. The same is true for free gluons. We evaluate the corresponding perfect quark-gluon vertex function, 
identifying in particular the "perfect clover term". First simulations for heavy quarks show that the mass is 
strongly renormalized, but again the renormalized theory agrees very well with continuum physics. Furthermore 
we describe the multigrid formulation for the non-perturbative perfect action and we present the concept of an 
exactly (quantum) perfect topological charge on the lattice. 



A large number of contributions to this con- 
ference are devoted to improved actions; there 
is no doubt that they are in fashion. This in- 
dicates a consensus that they represent a ray of 
hope for a great leap forward in lattice QCD. 
Most improvement procedures follow in one way 
or the other Symanzik's program [1], using per- 
turbation theory in the lattice spacing, e.g. [2]. 
Our work employs a different concept, utilizing 
renormalization group tools to construct an ap- 
proximately perfect action. A fixed point action 
(FPA) on a critical surface is an example of a 
perfect action, an action without any cutoff arti- 
facts. Using the fixed point action even at finite 
correlation length, a drastically improved scaling 
behavior has been observed for the 2d 0(3) model 
[3] and pure 4d SU(3) gauge theory [4], and is ex- 
pected also for other asymptotically free theories 
such as full QCD. 

1. Free fermions 

Fixed point actions have been constructed also 
for free fermions [5,6]. They can be obtained 
from iterating block renormalization group trans- 
formations (RGTs) with a finite blocking factor, 
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or more efficiently by sending the blocking factor 
to infinity and performing only one step. This 
amounts to a technique that we call "blocking 
from the continuum". It has been used exten- 
sively in a discussion of the Schwinger model [7] 
and for quarks and gluons [8]. One starts from 
a continuum theory, divides the coordinate space 
into lattice cells and defines lattice fields by in- 
tegrating over these cells. For free fermions the 
lattice action is then given by 

e -s[*,«] = J £ty£tyexp{- -0] (1) 

-\ E [*« - / K _ / *M dy ] }• 

Here i/>, ij) are continuum fields, s is the continuum 
action and c x is a unit hypercube with center x. 
S is the perfect action in terms of the lattice fields 
ty, ty. Finally a > is an arbitrary RGT parame- 
ter; for any choice of a all expectation values are 
invariant under this RGT. For a — > this is a 8 
function RGT, and a > "smears" the 8 function 
to a Gaussian. 

For fermions of mass m, this yields in momen- 
tum space the perfect action 

S[*,¥] = f d d p^!{-p)Af {p)-H{p) 

Jb 
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A'(p) 



h(p) 



n( P + 2tt/) 2 



Pl , = 2sin( Pl ,/2). (2) 



Here 5 =] — it, ir] d is the Brillouin zone, A-^ is the 
perfect lattice propagator and the function II is 
the Fourier transform of the step function which 
is 1 in Co and otherwise. This construction is 
illustrated in Fig. 1. The blocking attaches a II 
function to ^ and to ^P, and the points in the 
blocks are connected by the continuum propaga- 
tor. The sum over I probes all points inside the 
blocks. The II functions assure its convergence. 
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Figure 1. Construction of the perfect propagator 
for free fermions. 

For m = this is a fixed point action, and for 
general m we have an entire renormalized trajec- 
tory for free fermions. We can go back to coordi- 
nate space and write the perfect action as 

S[% *] = J2^^pA x - Whu + A(i - y)]*„.(3) 

It turns out that for any a > this action is local 
in the sense that p^ and A decay exponentially [5] . 
We can now tune a to optimize the locality. In 
d = 1 the summation in A-^ (p) can be done ana- 
lytically and the action has only nearest neighbor 
couplings iff 



a(m) 



(4) 



In higher dimensions we cannot achieve such "ul- 
tralocality" (i.e. a finite support for and A), 
but the above choice for a still provides an ex- 
tremely local perfect action, see Figs. 2 and 3. 
These figures show that we can also handle heavy 
quarks, where the action becomes even more lo- 
cal. 
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Figure 2. The exponential decay of |p^| on the 
4-space diagonal. 
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Figure 3. The exponential decay of |A| on the 
4-space diagonal. 

As an application to thermodynamics, we con- 
sider the following ratio for free massless fermions 



p 77T 2 

4 = = 0.3838..., 

T 4 180 



(continuum) 



(5) 



where p is the pressure and T the temperature. 
Fig. 4 shows this ratio at N t lattice points in the 4 
direction for Wilson fermions and the above fixed 
point fermion. For the latter the cutoff artifacts 
are strongly suppressed; the remaining artifacts 
at very small N t correspond to temperature de- 
pendent prefactors that we have neglected in the 
partition function when we preformed the RGT. 
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Figure 4. The ratio p/T 4 versus N t for free mass- 
less Wilson, fixed point and hypercube fermions. 



coupling 


m = 


m = 2 


Pi(1000) 


0.136846794 


0.0185415007 


Pi(HOO) 


0.032077284 


0.0031625467 


Pi(lHO) 


0.011058131 


0.0007898101 


Pl(llll) 


0.004748991 


0.0002501304 


A(0000) 


1.852720547 


0.8442376349 


A(1000) 


-0.060757866 


-0.0119736477 


A(1100) 


-0.030036032 


-0.0032647950 


A(1110) 


-0.015967620 


-0.0011445684 


A(llll) 


-0.008426812 


-0.0004622883 



Table 1 

Couplings for a free hypercube fermion action, 
which are perfect for 3-periodic configurations. 
Note that all coordinates of A and all non-/x co- 
ordinates of can be permuted and sign flipped 
arbitrarily, is odd in the fj, coordinate. 



To develop a practical scheme one must trun- 
cate the exponentially small tails in perfect ac- 
tions at a finite range. Here we consider a hy- 
percube action for quasi-perfect free fermions. In- 
stead of simply setting the couplings outside a 
unit hypercube to zero, we start from the obser- 
vation that one can easily obtain perfect actions 
in a finite volume with periodic boundary condi- 
tions by reducing the integral J g d d p in Eq. (2) to 
a discrete sum. In particular if we use 3-periodic 
boundary conditions, we obtain couplings for the 
free action that naturally live within a hypercube 
and the resultant hypercube Dirac operator is de- 
termined by nine couplings, which are given for 
two masses in Table 1. Using these couplings in 
the infinite volume represents an elegant trunca- 
tion scheme, which maintains the exact result in 
momentum space for p^ = 0,±27r/3. The ratio 
p/T 4 for this hypercube fermion (at m = 0) is 
also shown in Fig. 4. 

The perfect action reproduces the continuum 
spectrum exactly. Truncation to a hypercube dis- 
torts the spectrum a little. However, it is still very 
good. In Figs. 5 and 6 we plot the dispersion re- 
lation of the hypercube action for the massless 
and a massive case, respectively. We compare it 
to the standard Wilson fermion action, the D234 
action [9] (see also Fig. 4) and a recent proposal 
by the Fermilab group [10]. 

There is a major difference between the on- 



shell Symanzik improvement and the perfect ac- 
tion improvement after truncation. The on-shell 
improvement concentrates on obtaining the spec- 
trum correctly order by order in the lattice spac- 
ing (i.e., small momenta). On the other hand, in 
the truncated perfect actions, there are errors at 
each order of the lattice spacing. However, these 
errors are exponentially small in the range of the 
truncation so that the dispersion relation is good 
over nearly the entire Brillouin zone, as Fig. 5 
and 6 show. Still one may analyze specifically 
the errors for small momenta. For example in the 
case of massless free fermions the small momen- 
tum expansion of the Dirac operator is given by 

A'(p)- 1 = »7„p„[l + cip^ + c 2 p 2 ] + c 3 p 2 . (6) 

The breaking of rotational invariance depends on 
c\. In the Symanzik on-shell improvement, c\ is 
tuned to zero. For the Wilson fermion we ob- 
tain ci = —0.166 and for the hypercube action 
ci = 0.024. In the case of massive fermions the 
dispersion relation can be expanded in small mo- 
menta: E = m s + p 2 /2m,]. + where m s is re- 
ferred to as the static mass, and m^ as the kinetic 
mass. Ideally one has m s = m^ . In the Symanzik 
approach this can again be tuned. In the trun- 
cated perfect action m^ is close to m s but not 
exactly equal. In Fig. 7 we plot m^ as a function 
of m s for the hypercube action, and compare it 
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Figure 5. Comparison of dispersion relations at 
m = 0. 
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Figure 6. Comparison of dispersion relations at 
m = 2. 



with the Wilson action. The fact that is a 
little larger than m s for the hypercube fermion 
can also be seen in Fig. 6, where the curvature 
at p = is slightly too large. Tuning to the 
exact value improves the dispersion modestly at 
p << 1, but can make it much worse at p = O(l), 
as the Fermilab curve shows. 

The above FPA was obtained using an RGT 
which manifestly breaks chiral symmetry. Conse- 
quently, the lattice fermions resemble the Wilson 
fermions in their chiral properties. In fact, iterat- 
ing RGTs on the lattice, we obtain the above FPA 
if we start from Wilson fermions. However, per- 
fect actions with the properties of staggered fermi- 
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Figure 7. Static vs. kinetic mass. 

ons are also of interest. In particular the remnant 
U(1)®U(1) chiral symmetry simplifies the study 
of chiral symmetry breaking. In Ref. [11] it was 
shown how one can block staggered fermions from 
a fine to a coarse lattice (with blocking factor 
3) without mixing the pseudoflavors. Later this 
blocking prescription was used to construct a per- 
fect action for free staggered fermions [6,12]. The 
same results can be derived using the technique 
of "blocking from the continuum" [13]. 

The perfect action for free staggered fermions 
can be written in momentum space as 

S[X,X\ = j^yi f dtpXpi-pWrfW^XsW, 

Kp' (p) = ~ ia » (p) Kp' (p) + 0(p) 8 py ■ ( 7 ) 

The indices p, p' run over the sixteen pseudofla- 
vors. A unitary transformation connects this 
pseudoflavor basis to the Dirac ® flavor ba- 
sis. The matrix T^ p , (p) is unitarily equivalent 
to (7,, ® 1), and is given by ^(p)^.^ + 
Sp+fi^eWP-P")/ 2 , with the sign factor 7 7/1 (p) = 
(— 1)^"<J' P» . The functions a /1 (p) and (3{p) take 
the form 

™ =£ (pTWT^ n(p + 27rZ)2 + a - 

The RGT parameters c and a are used to optimize 



5 



locality of the action, a is a mass-like smearing 
parameter of the 8 function analogous to the pa- 
rameter a used before, and c is its kinetic coun- 
terpart. 

In the above representation the pseudoflavors 
are placed on a lattice with unit spacing. When 
gauge fields are introduced, it will be convenient 
to place the fermions on a finer lattice with spac- 
ing 1/2, distributing the pseudoflavors over the 
16 corners of the hypercube of the finer lattice. 
Then there is one variable per site of the finer 
lattice and the action takes the form 

s ix,x\ = ^2xx[iv/i{x - y)p/i{x - y) 

- y)\Xy (8) 

where p /1 (a; — y) is non-zero only when (x — y)^ £ 
{±1/2, ±3/2, . . .} and (x — y) v ^f, is an integer, 
and X(x — y) is non-zero only when (x — y) v is an 
integer for all v. Due to the symmetries involved, 
pi(x — y) and \{x — y) contain all the information. 

Again in d = 1 the summation over I can be 
done analytically and the action is ultralocal iff 

sinh(m) — m cosh(m/2) — 1 
a(m) = 1 , c(m) = . 

Following the example of Wilson fixed point 
fermions, we use the same parameters in higher 
dimensions for optimal locality. In Table 2 we 
give the largest couplings for mass m = and 
m = 2. 
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Figure 8. Dispersion relations for staggered fermi- 
ons at m = 0. 

Using the couplings for m = and omitting 



(x,y,z,t) 


pi(x + 0.5,y,z,t) 


X(x,y,z,t) 


m=0.0 






(0000) 


0.6617391 


0.0 


(0100) 


0.0441181 


0.0 


(0110) 


0.0046569 


0.0 


(0111) 


0.0004839 


0.0 


(1000) 


0.0234887 


0.0 


(1100) 


—0.0004933 


0.0 


(1110) 


—0.0009913 


0.0 


(1111) 


—0.0004819 


0.0 


(0200) 


0.0018423 


0.0 


(0210) 


0.0001419 


0.0 


^IzUU ) 


n nnnioi 1 
— U.UUUlzll 


u.u 


(1210) 


-0.0001011 


0.0 


m=2.0 






(0000) 


0.3586038 


0.9799873 


(0100) 


0.0154803 


0.0246899 


(0110) 


0.0012088 


-0.0015531 


(0111) 


0.0001282 


-0.0009711 


(1000) 


0.0071677 


0.0246899 


(1100) 


-0.0004502 


-0.0015531 


(1110) 


-0.0002537 


-0.0009711 


(1111) 


-0.0000832 


-0.0003371 


(0200) 


0.0003688 


0.0005581 


(0210) 


0.0000096 


-0.0000703 


(1200) 


-0.0000366 


-0.0000703 


(1210) 


-0.0000155 


-0.0000265 


Table 2 



The largest values of Pi(r) and A(r) for free stag- 
gered fermions at mass m = and m = 2. 

the distance 2 couplings, we obtain the dispersion 
relation for the free staggered fermions shown in 
Fig. 8. 

2. Free gauge fields 

Since gauge fields live naturally on the lattice 
links, the blocking scheme for free gauge fields 
is different from matter fields. An adequate pre- 
scription for their blocking from the continuum 
was given in [7,8]. One integrates over all straight 
connections of corresponding points in adjacent 
cells. Hence one identifies the lattice gauge field 
(living on link centers) 

a Ap) = ^ a "( p + 27rZ ) n "( p + 27rZ )(- 1 )''' 
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n„(p) = ^n(p) 



(9) 



where is a continuum gauge field. In analogy 
to Eq. (1) one finds an FPA for free gauge fields. 
In a general gauge it reads 

S[A " ] = J/pM-p^Up^Mp) 
» Ap > ~ » (p > V E P P?A?(p)-i 

Ag(p) = E ( f +27r02 J +« + lPl- (10) 

a and 7 are RGT parameters analogous to a and 
c in the case of staggered fermions. We tune 
both of them to optimize locality. The FPA for 
a = 1/6, 7 = —1/72 is in d = 2 the standard 
plaquette action, and in four dimensions it is still 
extremely local [8]. The perfect free gluon can be 
truncated to a unit hypercube in the same way it 
was done for fermions, and the dispersion relation 
for transverse gluons is still excellent, see Fig. 9. 
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Figure 9. Dispersion relation of the transverse 
gluon. 

To impose gauge invariance one can parame- 
terize this weak coupling FPA in terms of closed 
loops. In Sec 6. we suggest a 6 loop parameteri- 
zation of the gauge action (see Eq. 26), which is 
a good approximations for weak fields and which 
can therefore be used to initialize our multigrid 
minimization to obtain non-perturbative results. 



As we noted above the staggered fermions re- 
quire an RGT that distinguishes the pseudofla- 
vors. For compatibility, the blocking convolution 
of the gauge field also has to be altered. The 
reason is that in the interacting theory the gauge 
field has to couple to the pseudoflavors appropri- 
ately. This can be achieved by blocking the gauge 
field from the continuum in the following way 



J d'^M^a^x-y), (11) 



MAv) = { 







1 

M(y?) = { I- 2|y„ 




otherwise 

W\ < \ 
\ < \vA< ! 

otherwise 



Because of its shape we call the "Mansard" 
function. In momentum space the expression 
analogous to (9) reads 

AM = E a„(P+4irl)nf (p+4irl)(-l)'" 



K(p) = ^^n( P ). 

Pn 



(12) 



A fl (p) is 47r antiperiodic, because it lives on the 
link centers of a lattice with spacing 1/2. We 
just replace the function 11^ by 11^, which is the 
Fourier transform of the Mansard function, and 
the momentum p by p/2. The Brillouin zone B 
then refers to ]2w, 2w]. Also the fixed point prop- 
agator maintains the form given in Eq. (10), up 
to the same substitutions (and 2wl replaced by 
47rZ). Ultralocality in d = 2 and optimal locality 
in d = 4 now requires the momentum dependent 
smearing coefficients 

-oHGniwf). 



Note that a + jp^ is always positive, so that the 
RGT is well defined. 



3. Interacting theory 

If we introduce finite interactions, it is no 
longer possible to block the fields analytically. 
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One can either use perturbation theory or nu- 
merical RGTs. However, one may hope for some 
improvement in non-perturbative physics by us- 
ing just the (truncated) FPA of the free fermi- 
ons. A simple test is the mesonic dispersion re- 
lation after rendering the fermionic action gauge 
invariant by hand. In order to study the mesonic 
dispersion relation using the hypercube action, a 
quenched simulation with the Wilson gauge ac- 
tion was performed at (3 = 5.0 on a 6 3 x 18 lat- 
tice. We used the m = couplings of Table 1 
and we obtained a static meson mass of M = 3.0. 
Thus, there is a huge renormalization of the quark 
mass and the simulation is only relevant to heavy 
quark physics. Since these simulations did not 
use the perfect quark-gluon interactions this is 
acceptable. In heavy quark physics the kinetic 
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suggests that renormalization preserves Lorentz 
invariance in the truncated theory. 

Another important quantity for heavy quark 
physics is the Pauli term <? • B /2m,B- For a given 
lattice action, mj can be evaluated by introduc- 
ing a constant magnetic field [10]. The hyperfine 
splittings depend on this term. When the hy- 
percube action is made gauge invariant by hand, 
tub turns out to be clearly improved compared 
to the Wilson action, but it is still quite bad, see 
Fig. 11. In the Symanzik approach it is the clover 
term that is used to obtain the correct Pauli term. 
In the perfect actions this term arises naturally 
from the calculation of the perfect quark-gluon 
vertex. Making the free fermion action gauge in- 
variant by hand is only a crude approximation to 
a perfect action. Hence there is no reason to be 
surprised that we did not obtain exact results. 
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Figure 10. Meson dispersion relation at bare 
quark mass m = in a volume 6 3 x 18 at (3 = 5.0. 
The static meson mass M turns out to be 3.0, 
showing a large renormalization. 

mass of the meson plays an important role. If 
it is equal to the static mass M, one obtains for 
the effective speed of light \/ E 2 — M 2 /p = 1 for 
small momenta. In the continuum this condition 
is imposed by Lorentz invariance and it holds for 
arbitrary momenta. It can be used to probe the 
behavior of a lattice theory. In Fig. 10 we plot the 
result of simulations for heavy mesons. One sees 
a dramatic improvement for the hypercube fermi- 
ons compared to Wilson fermions. This result 



Figure 11. Static mass vs. "magnetic mass" tub- 
4. The perfect quark-gluon vertex 

We now discuss the interacting theory to first 
order in the quark gluon coupling g. Now the 
blocking from the continuum of a fermionic field 
involves the continuum gauge field. As an ansatz 
we write 

*«'(p) = ->P i {P+ 27rZ)n(p+ 2ttZ) (14) 

+ j£jdJ2 J d d qK fl (p+2wl,q + 2wl) X 
a;(p -q)\W(q + 2*l), 
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where A" are Hermitian generators of the SU(N) 
algebra. Gauge covariance requires the kernel 
function to obey 

[Pr-q^Kpfaq) = IL(p-q)IL(q) - IL(p). (15) 

In the framework of perturbation theory, it is 
most suitable to determine as follows: expand 
the continuum gauge phase in one cell in terms 
of dp , using boundary conditions that forbid any 
flux of dp across the surface of the cell. This 
expansion yields the "boundary condition kernel 
function" 

Kl c {v,q)= { -^^Y^%i[fiAv v -q v )fM, 

P f q f l v = \ 

fi„{Pv) = — TZTTo Vv 



Pi ~ 



(16) 



An alternative determination, which is more suit- 
able for the non-perturbative search of fixed 
points by a minimizer that performs inverse 
RGTs with block factor 2, proceeds as follows: 
block from the continuum to a lattice of spacing 
1 and also to a lattice of spacing 2. Relating the 
results by a block factor 2 RGT leads to a recur- 
sion relation and its iteration to the "recursion 
kernel" 



k„ [p, q) = ^2^ n(p/2") ^2^' n( 2^+r) 

:t)> ( 17 ) 



x sin 



^ o n(p/2") 
-jc ( - 

where 

K/i(P> «) = [ II cos ( Pv + wl »/ 2 ) x 

cos(g v -ir/ v /2)]/[l + r 2 ] , 

and I excludes the fj, component. In d = 1 the 
two kernels coincide, since both of them obey the 
condition (15). In d > 1 they are different, but 
extremely similar. Both have the small momen- 
tum expansion 



K,(P, q) = || [l + ^{p 2 ~ 4[p„(p„ - g„) + ql] 
— 5[p(p — q) + q 2 ]} + 0(momentum 4 )j 



where p = {Pn,p),q = {q,t,q)- They vanish at 
the origin, have a peak < 0.1 in the region where 
the momentum components are O(l), and drop 
quickly to zero at larger momenta. It turns out 
that the difference of the kernels is really tiny 
(O(10 -5 )) everywhere. 

If we apply our method of blocking from the 
continuum consistently to 0(gA fl ) we obtain the 
lattice action 

S[¥, ¥, Ay\ = S[¥, ¥] + S[AJ + V[9, ¥, (18) 
with an interaction term of the form 

v[% % A] = — / d d p d d q ¥(-p) X 
Jb 2 

gV ll (p,q)A;(p-q)X^(q). 

In momentum space we identified the vertex func- 
tion analytically [8]: 

V ll (p,q) = Af(p)- 1 AL(p-q)- 1 
x IL v {p + 2wl - q - 2wn) 

^ (p + 2irl — q — 2irn) 2 { 1 



II(p+27rZ) U(q + 2irn) 



i(tf+2-Kjl) + m i(rf+2mf) + m 
U(q + 2irn) 



+K v {p+2wl,q + 2wn) 
+K v {q + 2wn,p+2wl) 



2mf) + m 
II(p+ 2wl) 
2wjt) + m 



(p^-q^A^p-q)- 1 

EpiPp^^Kip-Q)- 1 

^A'^-i-A'tp)- 1 ]. 



(19) 



This formula might look a little complicated, but 
it becomes quite transparent if we consider the 
"Feynman diagrams" for perfect lattice perturba- 
tion theory, illustrated in Fig. 12. They refer to 
the Landau gauge in the continuum and explain 
essentially the lengthy term involving n . In- 
serting the general lattice gauge propagator A^ 
and adding the last two lines takes the vertex 
function to a general gauge. 

The action (18) is gauge invariant to O(g) since 
Vf, (p, q) obeys the Ward identity 

(p^qMp,q) = A'(g)- 1 - A'(p)" 1 . (20) 
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Figure 12. Construction of the perfect quark- 
gluon vertex function by building blocks which 
arise from blocking out of the continuum. 

To 0(A fl ) this action is free of all lattice arti- 
facts of 0(a n ) and 0(ga n ), where a is the lattice 
spacing. Possible artifacts come in at 0(g 2 a n ). 

The construction of the 3-gluon vertex and of 
higher order terms is straightforward from the 
same building blocks. Higher order terms include 
higher and higher dimensional sums in momen- 
tum space. They are all UV regular, since we 
integrate from the beginning over the very short 
distances (which brings in the II functions). How- 
ever, it becomes very hard to evaluate such terms 
and transform them to coordinate space. For the 
quark-gluon vertex discussed here - with its 2d- 
dimensional sum - this is still feasible in d = 4. 

5. The perfect clover term in coordinate 
space 

Let us first map the system on d = 2. To keep 
the number of parameters small and still observe 
locality as a crucial property, we even specialize at 
first on the case where the fermion fields are con- 
stant in the 2-direction. Then the inverse Fourier 
transform of the perfect vertex function can be 
performed analytically: 

X 

-**(! + *ff^M+i/ 2 )**+i 
-V x+l {l-igA lx+ll2 )y x 

+*^i(l + *ff^ M+i/2 )* x+i 



-*x+i« r i( 1 - < ff^M+i/ 2 )*4 

+g ^2 |*xO"2X(a:i - 2/1,2/1 - z x )iA 2 ^ z 

xyz 

+Vcr 3 0{x 1 -y 1 ,y 1 - Zl )F y y z }], (21) 

F being the gauge field strength defined on the 
plaquette centers. The interaction terms oc l,<r± 
are obvious by gauge covariance. The interesting 
parts are the functions x an d 0. They turn out 
to be ultralocal, which is also promising for the 
locality in higher dimensions. The complete table 
of their non vanishing elements at m = is given 
in [8]. 

m = 

_19_ _±_ 7 1.5 

120 120 120 120 




cm 



4.455 0.165 1.771 0.344 

m = 1 [icr 2 ] 

Figure 13. The 2d perfect clover term for fermion 
fields which are constant in one direction. The 
symbols x and O mark the positions of ^ and ^, 
and we show the plaquette couplings for m = 
(above) and m = 1 (below, in units [10 -2 ]). 

In particular 6 determines the "perfect clover 
term" in this case. We now concentrate on that 
term oc 03. An illustration of its plaquette cou- 
plings at m = and 1 is given in Fig. 13. We see 
that the usual clover term, given by the closest 
plaquette coupling in the case where ^, ^ are at 
the same position, is indeed dominant. However, 
there is at least one more contribution which is 
certainly not negligible. In the case where ty, 
are separated by one link, the coupling to the pla- 
quette attached to this link is suppressed by much 
less that an order of magnitude. A parameteriza- 
tion in terms of parallel transporters relates this 
term to the "staple". 

Now we come to the general 2d case and im- 
pose 3-periodic boundary conditions. Then the 
plaquette couplings do not follow uniquely from 
the link couplings - given in formula (19) - any 
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more. If we assume a behavior under unfolding 
as in the above, effectively Id case, we obtain the 
result described in Fig. 14 for m = 0.001 and 
m = 1. Again we observe that the "staple term" 
is only suppressed by about a factor 2 with re- 
spect to the usual clover term, whereas all the 
rest is much smaller. 



= 0.001 

[io- 2 ] 

3.742 
0.562 
0.205 



2.383 
0.755 
0.160 
0.008 









-i 













[io- 3 ] 

11.543 
1.119 
0.353 















> 


< — ( 


5 — 



6.701 
1.736 
0.353 
0.000 



plausible from the identity 



which displays this summed Pauli term as a nat- 
ural partner of the Laplacian. In general, how- 
ever, the summed Pauli term is a function of the 
fermion mass, namely 

s(m)= p) 2 [l-ll, (A = e ™-1). (22) 

This constraint is obeyed by the 2d plaquette cou- 
plings given in Figs. 13,14, and the summation 
of the general 2d case in only one direction also 
reproduces the result of the effectively Id case. 




0.461 Q 1.712 

-0.004 0.134 

-0.071 -0.060 

-0.057 * -0.091 




Figure 14. The plaquette couplings in the general 
2d perfect clover term after hypercubic trunca- 
tion. The darkest equivalence class has the largest 
coupling etc. (at m = 1). Again x and O mark 
the positions of ty, ty. 



If we sum up all these Pauli terms, i.e. if we 
consider the case where ty, are constant every- 
where, then we obtain | at m = 0. This becomes 




i 

Figure 15. The largest link couplings oc 7172 in 
the 4d perfect vertex function V\ after hypercubic 
truncation (in units [10 -4 ]). Also here x, O mark 
the positions of ^> , ^ . 

Finally we look at four dimensions. The largest 
couplings at m = 1 are represented in Fig. 15, 



11 



where we consider the clover terms oc 7172 in the 
vertex function V\. Contributions oc 7273 etc. 
are much smaller, and the rest follows by symme- 
tries. In all the cases where ^ , are separated by 
diagonals, the couplings are strongly suppressed. 
If they are separated only by one link, however, 
we observe in Fig. 15 the presence of important 
contributions, related to various types of staples. 
This confirms the trend that such staple terms are 
only suppressed by a factor ~ 2 with respect to 
the clover leaf. We therefore recommend includ- 
ing at least these extra terms in improved actions. 

In spite of the optimization of locality, we have 
probably still too many couplings for practical ap- 
plications. It turns out that a crude truncation 
in coordinate space destroys qualitative proper- 
ties such as the relation to lower dimensions by 
setting the remaining momenta equal to 0. So we 
do not choose that truncation scheme, although 
its results are quantitatively similar to the ones 
we show here. The small contributions should not 
just be omitted but projected onto the surviving 
couplings. This is exactly what periodicity did 
for us, but it's capacity is exhausted: a period 
smaller than 3 annihilates the entire kinetic term 
of the fixed point fermion and is therefore not 
sensible. 

A simple way to analyze the higher dimen- 
sional operators that are introduced by the per- 
fect quark-gluon vertex is an expansion in small 
momenta. Here we consider the terms in the ver- 
tex that have the structure 0",j(p,- — q{)Aj, which 
are the clover like terms. If we assume to be 
constant in Euclidean time (p4 =94), the leading 
term in small spatial momenta has the form, 

Vj{p, q,P4) = <rij{Pi ~ qi)[ij4{ci sinp 4 

+C2 sin 2p4) + do + di cosp4 — C2 cos 2p4]. (23) 

Remarkably, there are no higher order Fourier 
terms. The coefficients at various masses are 
given in Table 3. Thus, the perfect vertex sug- 
gests a more complicated form for the clover term 
than is usually used in the lowest order Symanzik 
improvement program. We already discussed the 
mass parameter mj for the hypercube action 
which was made gauge invariant by hand. If 
we include the information of the clover terms 
through couplings c±, do and d±, we improve 



m 


Cl 


C2 


do 


di 


0.01 
1 
2 


-0.219 
-0.0534 
-0.0107 


-0.0137 
-0.00287 
-0.000425 


0.206 
0.0676 
0.0183 


0.274 
0.711 
0.0150 



Table 3 

Coefficients of the clover like terms (described in 
the text) obtained from the small momentum ex- 
pansion of the quark-gluon vertex. 



tub considerably. This improvement is shown in 
Fig. 11. 

6. Non-perturbative perfect action 

So far we have emphasized our perturbative 
methods and results. The reasons are two fold: 
(i) Perturbation theory allows us to get analytical 
results, whereas the fully non-perturbative RGT 
with few exceptions (see next Section) must ul- 
timately be formulated numerically as a multi- 
grid minimizer. (ii) Any practical RGT scheme 
requires truncation in the range of the lattice ac- 
tion so the perturbative results are very useful for 
assessing the impact of this truncation and giving 
us a first approximation to the fixed point action. 
Nonetheless a fully non-perturbative perfect lat- 
tice action is the ultimate goal. Here we outline 
the basic steps and report on the current state of 
our research. 

The (quantum) fixed point action for lattice 
QCD, SlUf,, ty, can in principle be defined by 
integrating out continuum fields a^, ip, in blocks 
associated with fixed lattice fields, 

e -S[U,V,V] = J DcbiiD ^ D ^ App 6(G gf ) 

xe s[a ll , -0, V>] - TglU^a,,] - T f [¥, *, i>, -0] ^ 

where App is the Faddeev-Popov term and 
S{Ggf) fixes the gauge. As before, the gauge 
and fermionic blocking transformation terms, T g 
and Tf , define the particular RGT by correlat- 
ing the lattice fields with a local average over the 
continuum fields. We assume that the fermionic 
blocking term has been chosen to be quadratic 
in the quark fields. For the gauge blocking, we 
must project non-compact continuum fields onto 
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a compact lattice fields, but - as we will report 
in a separate publication - this requires non- 
Gaussian terms. The "classically perfect action", 
which corresponds to the limit h — > 0, is found by 
taking the minimum of the exponent on the right- 
hand side. The classical limit (or tree approxima- 
tion for the effective action) greatly simplifies the 
problem by removing the Faddeev-Popov ghosts 
and decoupling the quarks from the fixed point 
gauge action. Thus we are left with an extremum 
condition for the continuum gauge fields, 



S g [U ll ] = Min a ^[s g [a ll } + T g [U ll 



(25) 



J gV w i>.\ — L "gl^iH t J-gi^/i, "-in j • 

In general we still find it convenient to define 
the RGT in a fixed gauge internal to each block, 
so the minimization is done constrained to this 
gauge surface. Subsequently, the fixed point 
fermion action is found from the minimum of the 
quadratic form 



S/ [^,y>] + T / [iM,y;,y>] 



which can be solved by matrix inversion in the 
fixed background of the continuum gauge fields 
dp determined by Eq. 25. 

Our approach to this problem is to break it up 
into two steps. Use perturbation theory to ap- 
proximate the integrals from the continuum onto 
a fine lattice followed by a few iterations of the 
full non-linear RGT on a multigrid, as illustrated 
in Fig. 16. Let us describe this procedure in some 
detail for the gauge action. In the classical limit, 
the perturbative expansion of the perfect action 
is given as a sum over tree diagrams, with the 
0(g n ) term contributing uniquely to the 0(n — 2) 
monomial in the fields. So far we have computed 
only the lowest order quadratic contribution - the 
gluon propagator. From this we can start to ap- 
proximate the compact lattice action for smooth 
fields. For example, we have chosen a parameter- 
ization with 6 terms, (see Sec. 2) 



S[U] = 0.497 Q - 0.0268 



+ 0.0461 







+ 0.00544 ^jT^ 



+ 0.000432 



+ 0.00292 



which accurately reproduces the quadratic term 
for loops up to length 8. Each term is represented 
by a Manton action: e.g. Tr[Log(U)Log(U^)] 
where U is the appropriate Wilson loop for each 
term in the action (26). When we have deter- 
mined the 3-gluon vertex, a more accurate small 
field action will be formulated. 




blocking from 
continuum 



YJ/y/ ////////// 




(26) 



Figure 16. Non-perturbative perfect action by 
combined continuum and multigrid minimization. 

Now we re-express the RGT in Eq. (25) as a 
recursive blocking by factors of 2. Specifically, 
we choose to define our blocking transformation 
in the Landau gauge, expressed in a Manton like 
form. Each gauge link V f,(x) on the doubled lat- 
tice is related to the 16 pairs of product links 
u fl (y)u fl (y + p,) on the original finer lattice that 
connect the two hypercubes containing x and 
x + 2/x. The gauge blocking function is given by 

T g [Up,up] = Y J Y. Tr ^ H l( x 'y) H ^y)\ ( 27 ) 

x,/j y 

in terms of the Hermitian generators, 

H„{x;y) = -»£off[ETt(i)u M (y)u M (y + /*)], (28) 
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for each of the 16 double links. We have written 
and tested a code for an efficient minimizer for 
this problem. The multigrid minimizer enforces 
the gauge constraint, 

G x = J2 TrlLogiu^LogiuKy))] = 0, (29) 

inside each 2 4 block characterized by the point x, 
using Lagrange multipliers X x . By applying gra- 
dient descent in u fl (y) combined with a gradient 
assent in X r to 



■T[U, 



we obtain an efficient algorithm. Fortunately the 
gauge constraint does not present any essential 
difficulty to the formulation or application of the 
algorithm. Also the Manton form of the blocking 
term (27), since it is zero at its minimum, has the 
advantage that there is no need to subtract a field 
dependent factor to normalize the FPA. Our hope 
is that by an accurate perturbative starting point, 
the fixed point can be found with one or two ad- 
ditional iteration steps. Moreover it is also pos- 
sible within this general framework to simulate 
RGTs at finite correlation length by the MCRG 
to obtain quantum corrections to the classically 
perfect trajectory. In the next section we turn 
to a simple example where classical and quantum 
perfection actions can be computed analytically 
and compared. 

7. Quantum perfect vs. classically perfect 
topological charge 

In this Section we work out the concept of 
a quantum perfect (really perfect) topological 
charge and study the artifacts of the classically 
perfect charge [13]. As an example we consider 
the Id XY model, which is equivalent to a quan- 
tum mechanical particle on a circle. Its propaga- 
tor from angle at Euclidean time to (f> at time 
t can be written as a path integral 



-SeMr)] 



(<M|o,o) e = / d v {t) 

./(o,oH(<m) 

= J2 D<p(T)exp{-S[<p(T)]+iQe} 

n^TZ J (0,0)->(<b,t) 



where r £ [0,t], uq is the set of paths with the 
given endpoints and winding number Q, and we 
have included a 6 term. For fixed uq we can 
insert the free propagator on IR and obtain the 
Villain form 



(31) 




exp{-^(<£ + 27rQ) 2 +iQ0}. 



where J is the moment of inertia. Of course 
we can insert discrete points between and t 
and integrate over all possible intermediate states 
(Chapman-Kolmogoroff equation). For example, 
if we consider N equidistant points and insert the 
exact transfer matrix, then we obtain the perfect 
lattice description, 



<^o,o> e =(— ) Ens/., 



n jv j = 1 nj 



+ (01- <j>0 + 27TO!) 2 ] + iQO}, 



d(f>j 



(32) 



where Q = n\ + 7x2 + . . . + tin and At = t/N . 
nj is the winding number between the positions 
4>j- \ and The perfect action shown here is an 
RGT fixed point, since the Id FPA of a free scalar 
particle can be made ultralocal by suitable para- 
meters. However, this description not only con- 
sists of a perfect action but also of a prescription 
for the measure. We do not attach one particu- 
lar charge to a given lattice configuration, as it is 
usually done by some smooth interpolation. In- 
stead we integrate over all possible interpolations 
of the lattice configuration, including all winding 
numbers nj G ^ in each discrete step. This is 
the very idea of path integration. Now the topo- 
logical charge is given by a probability distribu- 
tions p{{4>j}, Q) for each lattice field configura- 
tion {4>j}- Thus it can never happen that "small 
windings fall through the lattice meshes". 

This is in contrast to the classically perfect ac- 
tion, where some topological charge gets lost in 
this way. This was observed to be significant on 
very coarse lattices in the 2d 0(3) model, CP(3) 
model and pure SU(2) gauge theory [14]. In 
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our case the classically perfect charge corresponds 
to the geometrical definition of the topological 
charge, i.e., neighboring lattice points are inter- 
polated by the shortest arc. Hence nj is confined 
to {—1,0,1}. Fig. 17 illustrates the artifacts of 
the classically perfect action. 

The energy eigenvalues of the Id rotor are 

E *w = h{ i+ hY> * = m,2,... 

and the correlation time is identified as 

* ■ . ... 



£i(0) - E o {0) 
The perfect topological susceptibility is given by 

x = \{Q 2 ) = -\d 2 e z t {e), 

where Z t {6) = (0, t\0, 0)# is the partition function. 
Perfect scaling requires the product x£ to be con- 
stant. Fig. 17 shows the scaling artifacts in this 
quantity for the standard and for the classically 
perfect action. 



0.053 



0.052 



0.051 



0.050 



0.049 



quantum perfect 

classically perfect 

standard 



1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 

Figure 17. Artifacts in the scaling quantity x£ as 
function of the correlation time £ in lattice units. 

The Villain type action is generally encoun- 
tered when one considers exact projections from 
continuum fields onto compact lattice fields along 
the lines discussed in Sec. 6. 

To summarize, we have outlined our program 
for constructing perfect actions for lattice QCD in 



the classical approximation. The results obtained 
from perfect lattice perturbation theory are suit- 
able for applications to heavy quark physics. 
Moreover they provide a basis for a multigrid 
construction of non-perturbative perfect actions, 
suited to strongly fluctuating fields. 
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